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Abstract 

Evolutionary dynamics in an uncorrelated rugged fitness landscape 
is studied in the framework of Eigen's molecular quasispecies model. 
We consider the case of strong selection, which is analogous to the 
zero temperature limit in the equivalent problem of directed polymers 
in random media. In this limit the population is always localized at a 
single temporary master sequence cr*(t), and we study the statistical 
properties of the evolutionary trajectory which cr*(t) traces out in se- 
quence space. Numerical results for binary sequences of length = 10 
and exponential and uniform fitness distributions are presented. Evo- 
lution proceeds by intermittent jumps between local fitness maxima, 
where high lying maxima are visited more frequently by the trajec- 
tories. The probability distribution for the total time T required to 
reach the global maximum shows a r~^-tail, which is argued to be 
universal and to derive from near-degenerate fitness maxima. The to- 
tal number of jumps along any given trajectory is always small, much 
smaller than predicted by the statistics of records for random long- 
ranged evolutionary jumps. 



* Permanent address. 
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"The concept of quasispecies is not just a model that involves any odd assumption; 
it shows how to view the darwinian world of replicating and mutating species from 
a physical viewpoint. " M. Eigen [j^ 

1 Introduction and motivation 

Eigen's quasispecies theory of molecular evolution is the simplest mathe- 
matical model that incorporates the central Darwinian paradigm of natural 
selection acting on variability created by random mutations. The model was 
originally developed to understand the conditions for the maintenance of in- 
formation in systems of self- instructive replicating macromolecules . Such 
systems can be realized in the laboratory in the form of populations of RNA 
strands which replicate in vitro in the presence of RNA replicase, display- 
ing a wide range of evolutionary phenomena P, ^, ^, The notion of a 
quasispecies refers to the structure of self-replicating populations, which 
typically consist of a distribution of related mutants centered around a most 
abundant master genotype (see below). The quasispecies structure plays an 
important role in the evolution of RNA viruses, where the presence of a wide 
range of mutants allows the virus to adapt rapidly to environmental changes 
P, ^ . On the other hand, the existence of an error threshold beyond which no 
localized quasispecies can be maintained (see Eq.(§)) places an upper bound 
on the genome length of RNA viruses [|10] . 

The mathematical structure of the quasispecies model has made it a fa- 
vored entrance way for statistical physicists into the field of biological evo- 
lutionQ. It was first observed by Leuthausser [|14| that the discrete time dy- 



namics (^ can be interpreted as a transfer matrix of a two-dimensional Ising 
model, where the genotype sequences become one-dimensional spin configu- 
rations that are coupled in the time direction through the mutation matrix 
([1|) |T5[. A similar relation can be established between (^) and the trans- 



fer matrix of a polymer directed along the time axis |]I6| In addition, 
Baake and coworkers have recently exploited the equivalence between quan- 
tum spin chains and a class of kinetic evolution equations closely related to 
the quasispecies model, in which mutation and selection occur in parallel 



12, 18 



In its most basic version, the quasispecies model is formulated in terms 



^Recent articles which review this connection are [^T], ^ |T^. 
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of standard chemical reaction equation^ written for the concentrations n„{t) 
of sequences a = (cti, ctat), each of which is composed of symbols drawn 
from an alphabet of K letters; the usual choice is a binary alphabet [K = 2), 
so that cTj takes the values and 1. The resulting sequence space consists 
then of 5* = 2^ points arranged on the vertices of an A^- dimensional hyper- 
cube. Each sequence a reproduces at a rate W{a), which may be taken as 
a measure of its fitness In the reproduction process errors occur with 
a mutation probability /i per site. The probability of creating a sequence a' 
when attempting to copy sequence a is therefore equal to 

Q^^(a a') = - nf-dnM (l) 

where 

N 

dH{<y,<y') = Y.(a,-a[f (2) 

is the Hamming distance between the two sequences, i.e. the number of 
digits in which the two differ. The dynamical evolution in discrete time is 
then given by 

n.(t+ 1) =Y,W{a')Q^{a' ^ a)n^,{t). (3) 

cr' 

Introducing the constraint of a fixed number of molecules leads to nonlinear 
loss terms on the right hand side of (^) . However since these can generally 
be transformed away, we ignore this complication here, at the expense of 
dealing with exponentially growing population numbers. 

The linear form of the evolution equation @ makes it plain that, for long 
times, the concentrations will approach that eigenvector of the evolution ma- 
trix W {a)Q ^{a a') which corresponds to the largest eigenvalue. Provided 
this eigenvector is localized in sequence space, it defines the quasispecies: A 
distribution of related mutants centered around the master sequence, which 
usually is the sequence with the maximum replication rate W{a). The most 
celebrated property of the model is its prediction of a sharp error threshold, 
where the quasispecies delocalizes, and the population spreads uniformly 
over sequence space. In terms of the sequence length and the mutation 
probability yU, the condition for a localized quasispecies takes the form 0, |^ 

< AT^ax = (4) 



^Similar equations arise in classical population genetics 
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where A denotes the selective advantage, a measure for the superiority of 
the master sequence compared to the other sequences. For the simplest case 
of a single peak fitness landscape, where the master sequence replicates at 
rate Wq and all other sequences replicate at rate Wi < Wq, the selective 
advantage is A = VFq/W^i, while for randomly distributed replication rates 



it is a functional of the rate distribution [|T^, 20|. In terms of the physical 



analogies described above, the error threshold phenomenon is equivalent to 



the thermal phase transition in the Ising model [jT4|, |T^, 0, ^ ^ and to the 
thermal unbinding of a directed polymer bound to an attractive columnar 
defect along the time direction |]16|, p!7| . 

Much less appears to be known about the evolutionary dynamics of the 
model, that is, the approach to the final quasispecies distribution from an 
initial localized or delocalized state. It was first pointed out by McCaskill 

that this dynamics should take the form of a "slowing optimization walk" 
through a succession of metastable states which correspond, in some sense, to 
local maxima in the fitness landscape. The separation of time scales between 
the (long) residence time in a metastable state and the (brief) transition time 
to the next maximum implies a punctuated pattern of evolution [^, |2^, |25| , 
which can be analyzed in analogy to variable range hopping in condensed 
matter physics p6|, ^ p^j . 

We should concede from the outset that the deterministic rate equations 
d^) are not an entirely appropriate description for this evolutionary regime, 
since the transition between two local maxima proceeds through the tails 
of the localized, metastable quasispecies, where the population numbers are 
small and fluctuations due to the finite number of molecules cannot be ig- 
nored 1^21, ^ It seemed nevertheless worthwhile to explore these ques- 



tions within the most basic, deterministic model, before turning to more 
sophisticated approaches. 

The present paper reports on some preliminary results from such an in- 
vestigation. To avoid the complications due to a finite error threshold, we 
consider a strong selection limit (to be described in the next section), in 
which the population is localized at a single site in sequence space at all 
times. This allows a direct comparison with simple schematic models of 
evolutionary dynamics, such as adaptive walks [^, Q and record 

dynamics |3^, ^ Adaptive walks describe the evolution of a genetically 
homogeneous population under the assumption that deleterious mutations 
(which decrease the fitness) are eliminated, while advantageous mutations 
spread instantaneously. The population then performs an uphill walk in 
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the fitness landscape, which terminates at a local maximum where no fit- 
ter one-step neighbors are available. In contrast, in the quasispecies model 
the population evolves through a chain of local fitness maxima which pro- 
gresses all the way to the global optimum. Some qualitative properties of 
these evolutionary trajectories are described in Section |, while Section | fo- 
cuses on a specific statistical feature, the total evolution time. A comparison 
with record dynamics is provided in Section |^, and some open questions are 
formulated in Section ^ 





Figure 1: Two evolutionary trajectories generated in two fitness landscapes 
with identical ordering of fitnesses but different fitness distributions. The 
monotonic increase in sequence number for the exponential distribution is 
fortuitous - there is no correspondence between the position of a sequence on 
the y-axis and its fitness 



Throughout the paper we will consider maximally rugged fitness land- 
scapes, in which the reproduction rate W{a) is chosen independently and 
randomly for each sequence. In contrast to simpler permutation invariant 
landscapes this makes it necessary to store all 2^ sequences during the 
iteration of the evolution equation^, restricting our numerical treatment to 

■^An approximate scheme which reduces the storage requirement from 2^ to N is de- 
scribed in |3ql. 
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rather short sequences; the results shown here are for = 10. A systematic 
analysis of the dependence on sequence length will be presented elsewhere 



m 



2 The strong selection limit 

The form of the strong selection limit is motivated by the analogy with 
the zero temperature limit in the associated problems of statistical physics. 
Following Peliti |jll], ^] we introduce an inverse "selective temperature" k > 
by writing the reproduction rates in the form 

Wia) = e'^^'^K (5) 

We want to take the strong selection limit A; — oo in such a way that only a 
few mutations occur in each time step. This requires to scale the mutation 
rate as 

/i = e"^^ (6) 

where 7 > is a constant. Inserting (^) and (P) into (^) it is clear that the 
sequence concentrations will grow for large k as 

n^{t) = e^^^'^'^l (7) 

In the limit k 00 the evolution equation @ then reduces to the recursion 

E{a, t + l)= max[E{a', t) + F{a') - -fdnia, a')]. (8) 

a' 

Since the term —jd^ suppresses mutations to far away sequences, we expect 
similar behavior for a model in which only nearest neighbor mutations are 
allowed, 

E{a,t + 1)= max [E{a' ,t) + F{a') - -fdHi(r,a')]. (9) 

dH(o-,o-')<l 

All results shown in this paper were obtained using the nearest neighbor rule 
(H), with the parameter 7 set to unityQ. 

We still need to specify the probability distribution p{F) of the fitnesses 
F{a). Two choices will be considered: The exponential distribution p{F) = 

, F > 0, and a uniform distribution on the interval [0, S], where S = 2^ . 
The reason for this particular scaling of the width of the uniform distribution 
will become clear below in Section ^ 

^For a discussion of the differences between the rules (§) and (||), see 
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3 Evolutionary trajectories 



It is evident from that, in the strong selection hmit k —* oo, the entire 
population resides at the global maximum of the function E{a,t). The po- 
sition of this maximum in sequence space will be referred to as the master 
sequence at time t, and denoted by cr*(t). At time t = the master se- 
quence is placed at a randomly chosen point u'-*^ by setting E{a^^\t) = and 
E(a,t) = —oo for a ^ a^^\ The subsequent time evolution a*{t) defines an 
evolutionary trajectory. 

Inspection shows that, after one or two time steps, such a trajectory 
passes exclusively through local fitness maxima, and eventually, after a total 
evolution time T, it invariably reaches the global fitness maximum. During 
the evolution, the master sequence spends increasingly long time intervals 
at local maxima of increasing fitness, with a few abrupt transitions in be- 
tween (Figure |I]). The number of transitions is small (see Figure much 
smaller than the number of local fitness maxima, which is on average equal 



to 2^/(A^ + 1) ^ 93 ig. This implies that most local maxima are bypassec 
by a typical trajectory. 
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100 

Number of visits at local maxima 



Figure 2: Fitness F{a) vs. the number of visits for all local maxima in a 
fitness landscape with exponential fitness distribution p{F). The particular 
landscape used here is near degenerate (gap size e ~ 0.02) 
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A quantification of tliis statement is shown in Figure ^, in which the fit- 
ness F of local maxima is plotted against the number of times it is visited 
by an evolutionary trajectory. These data were generated by going through 
all possible starting points a*^*^ in a fixed fitness landscape. The figure shows 
a roughly linear correlation between the fitness of a maximum and the loga- 
rithm of the number of visits. 



In relation to adaptive walks |30, 31, 32, 33|, which respond only to the 



relative ordering of fitnesses and not to their actual values, it is of inter- 
est to ask to what extent the set of maxima visited by a given trajectory 
is determined by the ordering of fitnesses. For this reason Figure |I| shows 
two trajectories evolving in landscapes which were generated using the same 
random numbers - thus having identical ordering of fitnesses - but with dif- 
ferent fitness distributions. It can be seen that the set of local maxima visited 
by the two trajectories is almost identical, apart from a small detour taken 
by the "uniform" trajectory, but the timing of the evolutionary transitions is 
markedly different in the two cases. With reference to G.G. Simpson's classic 
treatise [^, we may say that the fitness distribution affects only the tempo, 
but not the mode of quasispecies evolution. For a quantitative analysis of 
the temporal aspects we next turn to the distribution of evolution times. 



4 Distribution of evolution times 

Figure ^ shows the distribution P{T) of the number of time steps T required 
to reach the global fitness maximum, obtained by averaging over 500000 
landscapes with exponential and uniform fitness distributions. The time 
distribution for the exponential case displays a distinct maximum around 
T = 7, followed by a slowly decaying tail which is well described by the 
power law 

P{T) ~ T"^ (10) 

over roughly two decades. The distribution for the uniform case is much 
broader, but a similar power law tail can be seen for times T > 500. 

The power law (0) appears to be a simple consequence of the order 
statistics of uncorrelated fitness landscapes. Let F^^^ > F*^^) > ... > F^^^ be 
a realization of fitnesses arranged in decreasing order. As a measure of the 
spread in fitnesses among the most fit sequences we introduce the fitness gap 

e = - > 0, (11) 
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Figure 3: Distribution of the total time T required to reach the global fit- 
ness maximum. Symbols show data obtained by averaging over 500000 un- 
correlated fitness landscapes with exponential (+) and uniform (x) fitness 
distributions, while the full line shows the distribution ([T7|) obtained for 
record dynamics. The simulations with exponential fitness distribution were 
stopped after 1000 time steps 



which is a random variable characteristic of each fitness landscape. In the late 
stage of evolution the population will typically make a a transition from the 
second best sequence a^"^^ (or some local fitness maximum with comparable 
fitness^]) to the globally optimal sequence a^^^ . From the evolution rule (H) it 
is easy to see that, for small e, this transition will require a time of the order 
of 

T^a{N)/e, (12) 

where the coefficient a is determined by the early stages of the evolution 
Thus given the gap distribution Pg(e) the tail of the distribution 



38 



process 

of evolution times can be estimated to be 

P(T) ^ aT^^P^{a/T), 



(13) 



and a T ^ power law follows for T ^ a, provided that < Pg(0) < cxo. The 
gap distribution is given by -Pg(e) = e~'^ both for exponentially distributed 

'^In fact cr^^) is a local fitness maximum with high probability 1 — N/ (2^ — 1) w 0.990. 
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fitnesses, and for fitnesses distributed uniformly between and S, wlien S is 



large The striking difference between the two evolution time distribu- 
tions seen in Figure ^is related to the different scaling of the coefficient a{N) 
in (|12D with sequence length: In the exponential case a{N) ~ N^^'^, while in 



the uniform case a{N) ~ y/NS ||38| . 

To compute -Pg(O) for general fitness distributions, note first that the joint 
distribution of F^^^ and F^"^^ is given by ^ 



P2(F«, = S{S - (14) 

where Pc{F) = Jq dF' p{F') denotes the cumulative fitness distribution. The 
cumulative gap distribution is obtained by integration, 

Prob[F(^) - < e] = 
SiS - 1) / p(f(2))p,(f(2))^-2 / p(F(i)), (15) 

JO 

which tends to -Pg(0)e for e — 0. Thus we conclude that 

PM = S{S - 1) / dFp{F)'p,iFf-' (16) 



which is clearly finite and nonzero. 

The relationship (|l^) implies that the near-degenerate fitness landscapes, 
which have very small gaps, are the ones that give rise to anomalously long 
evolution times. Figure ^ illustrates this connection. The data shown as 
crosses were obtained from an average over exponential fitness landscapes, 
for which the fitness gap e was increased artificially by increasing the global 
fitness maximum according to F'^^^ F^^^ + 1. This is seen to immediately 
remove the power law tail (p^). 



5 Comparison to record dynamics 



A simple schematic analogue of the nonstationary (ever slowing) evolution 
process found in the quasispecies model is provided by the dynamics of 
records 



34, 35, 36, 39 



which is equivalent to evolution by long-ranged ran- 
dom mutations known in the classical literature as the theory of "hopeful 
monsters" [^, |30|. In the present context it reduces to the following rule 
for the motion of the master sequence o"*(t) in sequence space: At each time 
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Figure 4: Comparison of the distribution of evolution times for exponentially 
distributed fitnesses (+) with an ensemble of landscapes for which near- 
degeneracies (small gaps) have been removed (x). The latter data were 
averaged over 50000 realizations 



step, the population attempts a jump to another, randomly chosen sequence 
a' ^ a* it). The move is accepted, and a*{t + 1) = a', if F{a') > F{a*{t))] 
otherwise it is discarded and a*{t + 1) = cr*(t). Thus the current sequence 
cr*(t) represents the fitness record among the sequences which the population 
has encountered so far. 

Clearly this process gives rise to a step-like, punctuated pattern which is 
qualitatively similar to that shown above in Figure |l]. Here we are concerned 
with a quantitative comparison of statistical properties. Let us first com- 
pute the probability distribution of the total evolution time T for the record 
dynamics. Since the probability of finding the global fitness maximum in 
any jump is l/S, the probability that it has not been found up to time t is 
(1 — 1/Sy ~ e~*/'^ for large 5* and t. Taking the derivative one obtains 

P(r) = S'^e~^/^. (17) 

The typical evolution times are of the order of the number of sequences, much 
larger than in the quasipecies model. This demonstrates impressively the 
"guided" nature of quasispecies evolution [0, which is much more efficient 
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Figure 5: Probability distributions for the total number of evolutionary 
jumps required to reach the global fitness maximum. Symbols show data 
obtained by averaging over 500000 fitness landscapes with exponential (+) 
and uniform ( x ) fitness distributions, while the dotted line is the log-Poisson 
distribution ([T8[) predicted by record dynamics 



than a random search. Figure ^ shows how broad the distribution ([T7| ) is 
compared to that of the quasispecies model. Note, however, that for very 
long times (longer than S) Eq.(p!7D decays exponentially, faster than the 
degeneracy- induced power law (|T0|) . Taken literally, Eq.(|lOD implies that the 
mean evolution time is infinite. 

Next we consider the distribution P„ of the total number n of evolutionary 
jumps which occur on the way to the global fitness maximum. Adapting the 



results of Sibani and collaborators ||3g, for the case of record dynamics 



we find that P„ is a Poisson distribution with parameter In S, 

In Figure ^ this is compared to numerical data obtained for the quasispecies 
model. Again the distributions for the quasispecies dynamics are much nar- 
rower, showing that less transitions are required to reach the global max- 
imum. Simulations for longer sequences show that the mean number of 
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transitions increases sublinearly in A^, more slowly than the linear behav- 
ior predicted by (0) |8|, gg. 



6 Outlook 

The simplicity of the strong selection dynamics (||,^ suggests to use it for a 
dynamical characterization of different kinds of fitness landscapes. In con- 
trast to the random landscapes considered here, realistic fitness landscapes 
obtained e.g. from RNA folding contain extended neutral networks in se- 
quence space, in which the fitness (defined in terms of the RNA secondary 
structure) does not change . Central concepts of quasispecies theory have 



been extended to such landscapes Extended neutrality provides a dis- 



tinct mechanism for the appearance of punctuation patterns in evolution, 
since changes in the genotype do not show up in the phenotype, as long as 
the former moves within a neutral network [Q. 

Another interesting direction for further research inspired by the analogy 
with directed polymers is to include effects of environmental fluctuations, 
which amounts to making the fitness landscape time-dependent In the 



directed polymer analogy, the issue is the interplay between columnar disor- 
der, which is provided by the time-independent part of the landscape, and 



point disorder modeling the time- dependent variations |46|. It is well known 



that point disorder can depin a polymer from an attractive columnar defect 



in much the same way as thermal fluctuations This suggests the in- 

triguing possibility of an error threshold delocalization phenomenon induced 
by environmental fluctuations. 



Acknowledgements. I would like to thank T. Halpin-Healy and C. Karl 
for their contributions to this project, and K. Sneppen, H. Flyvbjerg and L. 
Peliti for useful discussions. This work has been supported in part by NATO 
within CRG.960662, and by NSF under Grant No. PHY99-07949. 

References 

[1] M. Eigen: Trends in Microbiology 4, 216 (1996) 
[2] M. Eigen: Naturwissenschaften 58, 465 (1971) 



13 



[3] D.R. Mills, R.L. Peterson, S. Spiegelman: Proc. Natl. Acad. Sci. USA 
58, 217 (1967). 

[4] CK. Biebricher: 'Darwinian Selection of Self-Replicating RNA 
Molecules'. In: Evolutionary Biology, Vol. 16, ed. by M.K. Hecht, B. 
Wallace, G.T. Prance (Plenum, New York 1983) pp. 1-52 

[5] O.K. Biebricher, W.C. Gardiner: Biophys. Chem. 66, 179 (1997) 

[6] J.S. McCaskill, G.J. Bauer: Proc. Natl. Acad. Sci. USA 90, 4191 (1993) 

[7] M. Eigen, J. McCaskill, P. Schuster: Adv. Chem. Phys. 75, 149 (1989) 

[8] E. Domingo: Clinical and Diagnostic Virology 10, 97 (1998) 

[9] R.V. Sole, R. Ferrer, I. Gonzalez- Garcia, J. Quer, E. Domingo: J. theor. 
Biol. 198, 47 (1999) 

[10] M. Eigen and C.K. Biebricher: 'Role of Genome Variation in Virus 
Evolution'. In: RNA Genetics. Vol. Ill: Variability of Virus Genomes, 
ed. by E. Domingo, J.J. Holland, P. Ahlquist (CRC Press, Boca Raton, 
PL, 1988) pp. 211-245 

[11] L. Peliti: 'Introduction to the statistical theory of Darwinian evolution'. 
|cond-mat/9712027 

[12] E. Baake, W. Gabriel: 'Biological evolution through mutation, selec- 
tion, and drift: An introductory review'. In: Annual Reviews of Gom- 
putational Physics VII, ed. by D. Stauffer (World Scientific, Singapore, 
2000) pp. 203-264 

[13] B. Drossel: 'Biological Evolution and Statistical Physics'. 
|cond-mat/0101409| (to appear in Adv. Phys.) 

[14] I. Leuthausser: J. Stat. Phys. 48, 343 (1987) 

[15] P. Tarazona: Phys. Rev. A 45, 6038 (1992) 

[16] S. Galluccio, R. Graber, Y.-C. Zhang: J. Phys. A 29, L249 (1996) 
[17] S. Galluccio: Phys. Rev. E 56, 4526 (1997) 

[18] E. Baake, M. Baake, H. Wagner: Phys. Rev. Lett. 78, 559 (1997) 



14 



[19] J.S. McCaskill: J. Chem. Phys. 80, 5194 (1984) 

[20] S. Franz, L. Peliti, M. Scllitto: J. Phys. A 26, L1195 (1993) 

[21] S. Franz, L. Peliti: J. Phys. A 30, 4481 (1997) 

[22] J.S. McCaskill: Biol. Cybern. 50, 63 (1984) 

[23] G.G. Simpson: Tempo and Mode in Evolution (Columbia University 
Press, New York, 1944) 

[24] CM. Newman, J.E. Cohen, C. Kipnis: Nature 315, 400 (1985) 

[25] S.J. Gould, N. Eldredge: Nature 366, 223 (1993) 

[26] W. Ebeling, A. Engel, B. Esser, R. Feistel: J. Stat. Phys. 37, 369 (1984) 

[27] J. Krug, T. Halpin-Healy: J. Phys. I France 3, 2179 (1993) 

[28] Y.C. Zhang: Phys. Rev. E 55, R3817 (1997) 

[29] D. Alvcs, J.F. Fontanari: Phys. Rev. E 57, 7008 (1998) 

[30] S. Kauffman, S. Levin: J. theor. Biol. 128, 11 (1987) 

[31] C.A. Macken, A.S. Perelson: Proc. Natl. Acad. Sci. USA 86, 6191 (1989) 

[32] C.A. Macken, P.S. Hagan, A.S. Perelson: SIAM J. Appl. Math. 51, 799 

(1991) 

[33] H. Flyvbjerg, B. Lautrup: Phys. Rev. A 46, 6714 (1992) 

[34] P. Sibani, M.R. Schmidt, P. Alstr0m: Phys. Rev. Lett. 75, 2055 (1995) 

[35] P. Sibani: Phys. Rev. Lett. 79, 1413 (1997) 

[36] P. Sibani, M. Brandt, P. Alstr0m: Int. J. Mod. Phys. B 12, 361 (1998) 
[37] P. Sibani, P.B. Littlewood: Phys. Rev. Lett. 71, 1482 (1993). 

[38] J. Krug and C. Karl (in preparation) 

[39] W. Feller, IntroducfAon to Probability Theory and Its Applications, Vol. 
2 (Wiley, New York 1971) 

15 



[40] H.A. David, Order Statistics (Wiley, New York 1970) 
[41] C. Karl: Diploma thesis (University of Essen, 2001) 
[42] P. Schuster: Biophys. Chem. 66, 75 (1997) 

[43] C. Reidys, C.V. Forst, P. Schuster: Bull. Math. Biol. 63, 57 (2001) 

[44] P. Schuster, W. Fontana: Physica D 133, 427 (1999) 

[45] M. Nilsson, N. Snoad: Phys. Rev. Lett. 84, 191 (2000) 

[46] I. Arsenin, T. Halpin-Healy, J. Krug: Phys. Rev. E 49, R3561 (1994) 

[47] L. Balents, M. Kardar: Phys. Rev. B 49, 13030 (1994) 



16 



